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PREFACE 

This  report  was  prepared  by  Dr.  Kazuhiko  Itagaki,  Research  Physicist,  Snow  and 
Ice  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  Funding  was  provided  by  DA  Project  4A161 102AT24,  Research  in 
Snow,  Ice  and  Frozen  Ground-,  Task  C,  Research  in  Terrain  and  Climatic  Con¬ 
straints  Environment-,  Technical  Effort  El,  Cold  Environment  Factors;  Work  Unit 
002,  Adhesion  and  Physics  of  Ice. 

Technical  review  was  provided  by  Stephen  Ackley,  Chief  of  the  Snow  and  Ice 
Branch,  and  Dr.  Motoi  Kumai,  CRREL. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  pur¬ 
poses.  Citation  of  brand  names  does  not  constitute  an  official  endorsement  or  ap¬ 
proval  of  the  use  of  such  commercial  products. 


ICING  RATE  ON  STATIONARY  STRUCTURES 
UNDER  MARINE  CONDITIONS 


Kazuhiko  Itagaki 


INTRODUCTION 

—'y  Icing  on  stationary  structures  such  as  oil  rigs  is  becoming  an  increasingly  serious  problem 
as  offshore  drilling  operations  in  the  subpolar  regions  become  more  common.  Vny  little,  if 
information  exists  on  this  subject.  Extensive  observations  have  been  made  of  icing  on 
the  upper  structures  of  moving  ships,  but  the  complexity  of  this  problem  makes  analysis  of 
the  results  very  difficult.  Even  the  generation  of  water  drops  in  this  case  involves  many  fac¬ 
tors,  such  as  windspeed,  wave  direction  relative  to  the  bearing  of  the  ship,  and  size  and  free¬ 
board  of  the  ship.  On  stationary  structures,  however,  the  problem  is  much  simpler,  since  the 
major  factor  in  drop  generation  is  whitecaps  produced  by  wind,  and  no  motion  of  the  struc¬ 
ture  is  involved. 

In  the  present  study,  a  theoretical  calculation  was  made  by  combining  the  data  available 
on  the  generation  of  drops  by  wind  with  data  on  the  proportion  of  ice  frozen  from  the  col¬ 
lected  water  .^Although  the  results  discussed  in  this  paper  are  preliminary  and  a  wider  range 
of  reliable  data  is  required,  the  general  trend  of  the  calculations  for  stationary  structures  is 
surprisingly  ^parallel  with  observations  made  on  board  ships. 


PROPORTION  OF  SEA  SPRAY  FROZEN  ON 
THE  STRUCTURE 

The  rate  of  ice  accumulation  R  can  be  calculated  by  the  following  formula: 

R  =  C,  C.  F  (1) 

where  Fis  mass  flux  of  the  water  drops  and  C,  and  C,  are  the  proportions  of  spray  frozen  on 
the  surface  (freezing  factor)  and  the  coefficient  of  capture  of  drops  respectively.  C,  can  be 
close  to  unity  for  larger  drops  such  as  sea  spray.  Although  many  other  factors  may  contri¬ 
bute,  Cf  seems  to  be  a  strong  function  of  the  air  temperature.  The  relationship  between  mean 
air  temperature  and  freezing  factor  (frozen  mass/captured  mass)  obtained  from  data  given 
by  Ono  (1964),  and  reproduced  in  Table  1,  indicates  that  C,  increases  linearly  as  the  air  tem¬ 
perature  decreases,  as  shown  in  Figure  1 . 


Table  I.  Amounts  and  cblorinities  of  accumulated  brine  (after  Ono  1964). 

Pareatbeses  ladkate  calculated  values  assuming  ckloriiity  of  respective  value  of  sea  spray  (19^*»). 
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Figure  1.  Freezing  factor  as  a  function  of  air  temperature  us¬ 
ing  data  of  Ono  (1964).  Numbers  next  to  the  data  points  in¬ 
dicate  the  amount  of  sea  spray  collected  by  the  icing  gauge 
(g/hr).  Solid  line  is  simpler  regression  equation. 
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The  linear  regression  analysis  of  the  air  temperature 
versus  C,  indicates  that  the  regression  line  intersects 
the  jr-axis  at  0.48  °C.  The  x-intercept  is  expected  to  be 
between  O.C°  and  the  freezing  point  of  seawater 
(-1.9°C)  unless  the  air  is  very  dry.  The  wide  scatter  in 
the  data  points  may  be  the  major  source  of  the  discrep¬ 
ancy  in  the  regression  line  intercept.  C,  may  also  be  a 
function  of  the  mass  flux  F,  making  the  data  points 
more  widely  scattered.  In  addition,  a  possible  nonline¬ 
arity  in  temperature  versus  C,  was  disregarded  for  the 
regression  analysis.  Since  data  points  were  scattered  so 
widely,  the  simpler  form 

(2) 

was  used  instead  of  the  more  complex  equation,  bringing  the  x-intercept  to  0.48  °C. 

This  regression  line  is  shown  as  the  solid  line  in  Figure  1 .  By  use  of  eq  2,  C,  becomes  unity 
at  -12.878  °C  (T,),  but  again  this  temperature  could  be  considerably  higher  than  the  actual 
case.  Brine  pockets  in  sea  ice  freeze  at  about  -18  °C,  limiting  the  lowest  temperature.  Brine 
drains  out  faster  from  the  accreted  ice  than  from  sea  ice  because  of  its  porous  structure; 
therefore,  T ,  for  brine  would  not  be  as  low  as  that  for  sea  ice.  More  accurate  data  are  needed 
to  clear  up  these  points.  The  calculated  C,  values  as  a  function  of  temperature  are  shown  in 
Table  2. 


MASS  FLUX 

Mass  flux  can  be  written  as 

F  =  j  n(r)  V  r>dr  =  q\  n(0)  V<t?d<t>  (3) 

where  n(r)  is  the  number  of  drops  of  radius  r  (diameter  <t>)  in  unit  volume,  V  is  the  wind 
velocity  and  q  is  the  density  of  water.  n(r)  is  a  function  of  wind  velocity  and  height  of  obser¬ 
vation.  In  icing  of  a  ship,  the  mass  flux  is  a  comglex  function  of  the  wave  direction  and 
height,  the  wind  direction  and  velocity,  the  ship’s  bearing  and  speed,  and  the  overall  struc¬ 
ture  and  size  of  the  ship.  In  icing  of  a  stationary  structure,  the  mass  flux  at  a  fixed  height  is 
primarily  dependent  upon  the  windspeed. 

Lai  and  Shemdin  (1974)  compiled  some  data  from  other  researchers  (Toba  1961,  Monahan 
1968)  together  with  their  own  data  on  the  size  distribution  of  sea  spray  in  their  Figure  12 
(reproduced  and  shown  in  Fig.  2  here).  Generally,  the  distribution  of  the  number  density  is 
inversely  proportional  to  the  square  of  the  drop  diameter  &  as 

n(<t>,Y)  =  (4) 

where  the  constants  A(V)  were  determined  for  each  windspeed  plot  for  a  diameter  of 
4>  =  10'1  cm  (listed  in  Table  3). 

These  constants  A(V)  were  adjusted  for  a  height  of  13  cm  from  the  average  water  surface 
for  the  extrapolated  windspeed  at  10  m  from  the  water  surface  in  the  original  Figure  of  Lai 
and  Shemdin  (1974).  The  mass  flux  F  was  then 


Table  2.  Temperature  ver¬ 
sus  freezing  factor. 
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temperature  (°C) 
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0.1 55 
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0.621 

-10.000 
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-12.000 

0.932 

-12.878 

1.000 

C,  =  0.07765  T 


-  % 


d(fim) 


Figure  2.  Size  distribution  of  sea  spray  [from  Lai  and  Shemdin, 
Journal  of  Geophysical  Research,  79(21):  3055-3063,  1974; 
copyrighted  by  the  American  Geophysical  Union]. 


F-  (5) 

Since  the  drop  diameters  observed  by  Lai  and  Shemdin 
(1974)  were  between  50  and  700  nm,  integration  in  this 
domain  yielded  F  as 

F=  qV  A(V)  [(7 x  10~4)J  - 

(5  x  10'5)*]  kg/m's.  (6) 


Table  3.  A(V)  and  windspeed  ver¬ 
sus  mass  flux,  calculated  from 
Lai  and  Shemdin  (1974). 


Wind 

(m/s) 

A(V)  at  I  mm 
(no./m't 

Mass  flux 
(kg/m's) 

12.6 

2.0  x  10-' 

3.21  x  10 

IS. 8 

4.5 

9.07 

16.4 

5.9 

12.35 

18.9 

24.0 

57.89 

The  experimental  data  points  of  Lai  and  Shemdin  and  of  Toba  seemed  systematically 
deviated  from  the  Held  data  points  of  Monahan  for  a  windspeed  that  is  5  m/s  lower.  There¬ 
fore,  all  the  experimental  data  points  for  windspeed  were  shifted  5  m/s  lower  to  maintain 
consistency  of  the  results. 
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Figure  3.  Proportional  factor  of  number  density  A(V)  as  a  function  of 
windspeed. 


A  third-degree  polynomial  regression  fitted  the  relationship  between  windspeed  V  and 
A(V)  fairly  well  (Fig.  3): 

A(V)  »  A<,  +  A>V+AlV2  +  A)V’  (7) 

where  the  constants  were  A0  =  -53.5173,  A,  =  11.3119,  A,  =  -0.7943  and  A,  =  0.01864. 
Mass  flux  F  was  calculated  and  is  shown  in  Table  3  for  the  corresponding  windspeed. 


ICE  ACCRETION  RATE 

The  ice  accretion  rates  R ,  calculated  from  the  above  data  using  eq  1,  are  shown  in  Table  4 
at  2  °C  intervals,  using  a  capture  coefficient  C«  of  1 . 

In  order  to  obtain  the  semi-empirical  relationship  between  temperature  T  and  windspeed  V 
for  a  fixed  ice  accumulation  rate  R,  eq  2,  6  and  7  were  combined  into  eq  1  and  solved  for  T: 

t  _ _ AtR _  /o\ 

VAo+  FM.+  VnA,+  FM, 

where  A<,-A,  are  the  same  as  in  eq  7  and  A,  =  12.8781.  The  results  are  shown  in  Figure  4 
with  R  as  a  parameter.  Since  the  flux  was  calculated  from  values  observed  near  the  sea  sur¬ 
face  (13  cm  above  the  surface),  direct  comparison  with  ship  icing  data  obtained  from  on¬ 
board  observations  is  difficult.  However,  the  general  trend  showed  surprisingly  good 
parallelism  with  the  diagrams  given  by  previous  authors. 


DISCUSSION 

The  icing  rates  measured  onboard  ships  by  Ono  (1964),  using  an  icing  gauge  (collecting 
cross  section  of  2.2  by  14.5  cm),  are  shown  by  dots  in  Figure  4.  The  rates  were  converted  into 
kg/m1  hr  and  are  shown  by  the  numbers  next  to  the  data  points.  It  is  difficult  to  compare 
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Table  4.  Calculated  ice  accretion  rate. 


Wind 

(m/s) 

Air 

temperature 

CO 

Ice 

accretion  rate  R 
(kg/m'hr) 

Wind 

(m/s) 

Air 

temperature 

CO 

Ice 

accretion  rate  R 
(kg/m’hr) 

12.6 

-  2.000 

0.180 

16.4 

-  2.000 

0.690 

0.360 

1.381 

0.539 

-  6.000 

2.071 

0.719 

2.762 

-10.000 

0.899 

-10.000 

3.452 

-12.000 

1.079 

4.143 

-12.878 

1.158 

-12.878 

4.446 

15.8 

0.507 

18.9 

-  2.000 

3.237 

1.015 

6.473 

1.522 

9.710 

2.029 

12.947 

2.537 

16.183 

3.044 

-12.878 

3.267 

-12.878 

20.841 

these  points  directly  with  the  calculated  values  because  the  major  mechanism  of  sea  spray 
generation  in  Ono’s  case  is  not  the  wind  but  rather  the  pounding  of  the  waves  against  the 
ship’s  hull.  Generally,  the  observed  icing  rates  found  by  Ono  correspond  to  the  calculated 
value  in  Figure  4  at  5  to  8  m/s  higher  windspeeds. 

Figure  5,  from  Tabata  (1966),  depicts  the  relationship  of  the  icing  class  to  the  temperature 
and  relative  windspeed  for  350-ton  and  450-ton  class  cutters.  It  is  clear  that  increasingly 
higher  windspeed  is  required  to  produce  considerable  icing  above  -5  °C,  while  little  tempera¬ 
ture  effect  is  seen  below  -5  °C.  The  same  trend  can  be  seen  in  Figure  4,  although  a  much 
higher  windspeed  is  needed  for  the  stationary  structure  to  produce  similar  icing  rates. 

Total  ice  accumulation  rate  was  measured  on  a  Japanese  cutter  by  Tabata  et  al.  (1963)  and 
correlated  with  air  temperature  and  Beaufort  wind  scale  as  shown  in  Figure  6.  The  area  of  ice 
coverage  was  roughly  150  mJ  and  the  maximum  accumulation  rate  at  Beaufort  wind  scale  5 
and  -7°C  was  13  kg/mJ  hr  (1.99  tons/hr),  which  corresponds  to  a  windspeed  of  19  m/s  or 
Beaufort  8  in  the  present  calculation.  The  maximum  accretion  rate  obtained  by  the  icing 
gauge  (Ono  1964)  was  61  kg/mJ  hr  for  -7°C  and  15.4  m/s,  which  is  again  at  least  one  order 
of  magnitude  larger  than  in  the  present  calculation.  The  difference  in  the  spray  production 
mechanism  on  ships  compared  with  that  on  stationary  structures  may  account  for  the  differ¬ 
ence  in  icing  rates. 

The  mass  flux  observed  under  laboratory  conditions  (Toba  1961,  Lai  and  Shemdin  1974) 
corresponded  to  that  under  natural  conditions  (Monahan  1968)  with  a  windspeed  of  5  m/s 
less.  Two  major  reasons  possibly  contributed  to  this  discrepancy.  First,  the  vast,  open  wind¬ 
ward  water  (fetch)  creates  roughness  considerably  higher  than  laboratory  conditions.  Sec¬ 
ond,  the  humidity  near  the  surface  is  near  saturation  in  natural  conditions,  whereas  a  less  hu¬ 
mid  environment  was  adopted  in  the  laboratory.  Therefore,  the  laboratory  data  were  ad¬ 
justed  5  m/s  lower  to  obtain  more  realistic  figures. 

Lai  and  Shemdin  (1974)  measured  the  height  distribution  of  the  drops  and  found  no  drops 
above  the  height  Z0<,  depending  on  the  windspeed.  However,  these  measurements  are  diffi¬ 
cult  to  apply  for  the  natural  case  because  of  the  differences  in  surface  roughness  and  humidi¬ 
ty.  Furthermore,  the  freezing  factor  may  be  different  for  the  higher  flying  drops  because 
evaporation  may  increase  salt  concentration  in  the  drops. 

The  effect  of  fetch  on  the  flux  droplets  and  Z0d  was  studied  by  Wang  and  Street  (1978). 
They  indicated  in  their  Figure  7  a  rather  parabolic  increase  of  Ztd  with  increase  in  fetch.  No 


Air  temperature.  °C 


Figure  5.  Effect  of  air  temperature  and  relative  windspeed  on  ice  accretion  class  (after 
Tabata  1969)  (copyright  T.  Tabata;  reprinted  by  permission). 
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Figure  6.  Total  ice  accumulation  rate  measured  on  the 
Japanese  cutters  Chitose  and  Yuburi  (after  Tabata  et  al. 
1963)  ( copyright  T.  Tabata;  reprinted  by  permission). 


saturation  was  attained  within  their  limit  of  fetch  (12.S  m).  Because  no  appreciable  differ¬ 
ence  in  droplet  size  distribution  can  be  observed  in  their  Figure  6,  which  shows  a  droplet  his¬ 
togram  for  different  fetches,  mass  flux  at  the  fixed  height  increases  as  the  fetch  increases.  It 
is  difficult  to  extrapolate  their  experimental  results  into  the  natural  condition.  However,  a 
5-m/s  difference  in  windspeed,  corresponding  tc  the  difference  in  mass  flux  between  natural 
observations  and  laboratory  experiments,  seems  reasonable. 

Neglected  factors  such  as  sea  spray  flux,  windspeed  and  water  temperature  may  be  re¬ 
sponsible  for  the  wide  scatter  of  the  data  points  in  Figure  1 .  Actually,  most  of  the  low  sea 
spray  flux  data,  as  indicated  in  this  figure,  lie  above  the  regression  line,  indicating  that  more 
factors  have  to  be  taken  into  account.  However,  limitations  in  the  available  data  precluded 
further  refinement. 

Considering  these  differences  in  icing  conditions,  especially  differences  in  the  drop  genera¬ 
tion  mechanisms,  direct  comparison  between  the  data  obtained  from  ship  observations  and 
from  the  theoretical  calculations  described  here  for  stationary  structures  may  be  inappropri¬ 
ate.  The  general  trend  found  in  Figures  S  and  6,  obtained  from  ship  observations  made  by 
Tabata  (1969)  and  Tabata  et  al.  (1963),  shows  strong  similarities  between  their  data  and  the 
theoretical  results  shown  in  Figure  4,  however.  Those  similarities  seem  to  indicate  that  the 
approach  of  the  present  analysis  is  in  the  correct  direction,  but  that  there  is  a  need  for  more 
accurate  icing  data  from  stationary  structures  to  compare  with  this  particular  analysis. 
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